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Abstract 

We present a new thermodynamical approach to multiparticle production in high 
energy hadronic interactions, making use of the formalism of infinitely divisible 
power series distributions. This approach allows us to define new observables, 
linked to the system fugacity, which characterise different classes of events. 
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1 Introduction 



The phenomenological analysis of many-particle final states in hadron-hadron collisions in the 
GeV region has been successfully carried out[|l]] using a two-component model: each event 
is assigned to one of two classes, called 'soft' and 'semi-hard', which correspond to events 
without mini-jets and events with mini-jets, respectively. We also assume that the multiplicity 
distribution (MD) in each class is described by a negative binomial (Pascal) distribution (NBD), 
of course with different parameters n, k in each class. This model was successful in describing 
the shoulder in MD's, the oscillations of high rank moments thereof^ and also forward back- 
ward multiplicity correlations.^] There are also experimental indications that the two classes 
behave differently in the TeV region. [|3p 

In order to extrapolate to the LHC region we have to make some assumptions on the be- 
haviour of the parameters of the two NBD's which we summarise as follows: 

* The overall average multiplicity grows as In 2 y/s. 

* The soft component average multiplicity grows as In y/s and the MD obeys KNO scaling 
(thus fc so ft is approximately constant.) 

* The average multiplicity in the semi-hard component is approximately twice as large as 
in the soft one. 

-k Three scenarios have been examined for the behaviour of k semi . haTd : 

1. same behaviour as for the soft component, i.e., it is constant (therefore KNO scaling 
is satisfied); 

2. fc~mi-hard ~ m V^> implying a strong violation of KNO scaling; 

3. it follows a QCD-inspired behaviour; KNO scaling is attained only asymptotically. 
This last scenario is intermediate between the first two. 

Of course, using NBD's means we can explore the clan structure: [|J |J the average number 
of clans N and the average number of particles per clan are defined by 

jv = *ta(i + -) i «.= fcln(1 + a/t) . 0) 

It turns out that the second and third scenarios show a number of clans which is rapidly 
decreasing with cm. energy (accompanied by a fast increase of the average number of particles 
per clan). This is surprising, and in this paper we try to understand the implications of this result 
at parton level using thermodynamical concepts. 

To connect the hadronic and partonic levels we use the generalised local parton-hadron 
duality (GLPHD),^] which says that all inclusive distributions are proportional at the two levels 
of investigation: 

Vn) = P* Qn,partons (yi,...,y n ), (2) 
which corresponds for NBMD parameters to 

^hadron ^partem; 'T'hadron P ^parton' (3) 

GLPHD will be applied separately to soft and semi-hard components. 
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2 A new thermodynamical approach 



The thermodynamical approach to multiparticle production has a long history which cannot be 
summarised here. We would just like to attract the reader's attention to the result[[7p that, to 
leading order in the allowed rapidity range, the generating function (GF) for the MD has the 
form of an infinitely divisible distribution (IDD). 

Keeping in mind this result, we propose the following approach. 

The partition function in the canonical ensemble, Q n (V, T), for a system with n particles, 
volume V and temperature T, is linked to the partition function in the grand-canonical ensemble 
with fugacity z, Q(z, V, T), by the well known relation 

Q(z,V,T) = J2z n Qn(V,T). (4) 

n 

Quite in general, in the grand-canonical treatment the probability of finding n particles in a 
system is given by 

v(n) = sfeW <5) 

That is to say, for a thermodynamical system the MD belongs to the class of power series 
distributions (PSD's), and is characterised indeed by the following form: 

t \ anb " (e.\ 
= ^ (6) 

with constants a n , b. 

We therefore propose, given a MD in power series form, the following correspondence with 
Eq. ©: 





< — 


7 

> 


I 





(7) 

Q< — >l(b) =p(0)- 1 . 

When the PSD is also IDD, then we know it can be cast in the form of a compound Poisson 
distribution, such that 

p(0) = e^. (8) 

In a two-step approach, N is average number of objects (clans) generated in the first step. 
This way of describing the partonic cascade is well known: []5p the ancestors (first step) are 
independent intermediate gluon sources; it is their thermodynamic properties which we want to 
explore. 

In our thermodynamical approach, N becomes of fundamental importance since Eq.s (0) 
and (||) imply 

N = -lnp(O) = InQ; (9) 

all thermodynamical properties can be obtained by differentiating N. 

From the standard relation PV = k B T\n. Q, we obtain the equation of state 

PV = Nk B T, (10) 

which says that clans form a classical ideal gas. 
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The negative binomial (Pascal) distribution belongs to both classes, power series and IDD. 
The standard form, from which the correspondence with the partition function can be obtained, 
is the following: 

. > k(k + l)...(k + n- 1) / n V / k \ k 

The identification we propose in our approach is 

k(k + l)...(k + n-l) 
tin j j 

• 5 (12) 



n + A;' 

7(6) = 



n + k 



Notice that fc(V, T) is the canonical partition function for a system with 1 particle; it is in 
our approach an unknown function of V and T. 
Finally notice that b is the fugacity z: 

* = ?rb' (13) 

When the ancestors are created early in the evolution, at larger virtualities and with higher 
temperature, they tend to follow a quasi-classical behaviour, as the production of a new ancestor 
is competitive with the increase in gluon population within each clan. This results in a relatively 
large value of the k parameter, i.e., a small amount of aggregation. When the number of partons 
per clan is very small (close to 1; k is very large) then essentially each parton is a clan, and the 
equation of state reduces basically to that of an ideal gas (quasi-classical behaviour): 

PV « nk B T. (14) 

Via GLPHD, we expect a similar situation to hold at hadron level. This behaviour is qualita- 
tively close to that of soft events as well as of scenario- 1 semi-hard events. 

When the ancestors are created later in the evolution, at lower virtualities and with lower 
temperature, they tend to remember their quantum nature, as newly produced gluons prefer to 
stay together with other clan members rather than initiate a new clan. This results in a relatively 
small value of the k parameter, i.e., a larger aggregation and larger two-particle correlations. 
When the number of partons per clan begins to grow, the equation of state for partons becomes 
more and more different (quasi-quantum behaviour), but that for clans remains that of an ideal 
gas. 

PV = Nk B T = kln{l + n/k)k B T. (15) 

Via GLPHD, at hadron level we recognise the behaviour of scenario-2 and scenario-3 semi-hard 
events. 

It is interesting now to calculate some thermodynamical quantities. 

The Helmholtz free energy can be rewritten in a form symmetric in n and k: 

A , / k 



knT 



fi/i- PV = n\n(^l + ^ +k\n(l + ^y (16) 
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The average internal energy is 



U , ^ 2 fdN\ - fdlnk 



The entropy is 



n 



which coincides with —A/T in the limit of (dk/dT) v — > 0, since also U — > 0. 
For further discussion of thermodynamical quantities, see Ref. [Rp. 



3 Clan behaviour as a function of fugacity 

Relying on GLPHD, we analyse first experimental data on the fugacity and the related a param- 
eter: the NBD satisfies the recurrence relation 



(n + l)p(n + 1) 



p(n) 

where 

nk 



a + 6ft, (19) 



(20) 



ft + k n + k 

From Eq. dl2| ) it is seen that b is the fugacity. In Figure [j] we show for each component and 
each scenario the energy variation of the parameters a and b. The points come from NB fits 
to experimental MD's, the lines show the predictions from the extrapolation mentioned in the 
introduction. 

The a parameter corresponds to the average multiplicity for a classical (Poisson) system. 
The relative behaviour of b and a = kb as the cm. energy increases can be considered an 
indication of the relative importance of a behaviour closer to a quantum one, i.e. harder, with 
respect to a behaviour closer to a quasi-classical, i.e. softer, for a class of events. A very slow 
increase of b with cm. energy and an almost constant behaviour of a is the main characteristic 
of the class of soft events and of scenario- 1 semi-hard events. A very fast decrease of a in 
scenarios 2 and 3 and larger values of the fugacity b characterise harder events: the assumption 
of strong KNO scaling violation for the semi-hard component (an extreme point of view with 
respect to that of scenario 1) implies a completely new panorama. 

Then we explore the dependence of clan parameters on the fugacity b, induced by its energy 
evolution: 

N = k\n(l + y) =-k\n(l-b); (21) 



k 

n = n = b - (22) 

c kln{l + n/k) (b - 1) m(l - b) ' V ; 

Notice that the average number of particles per clan only depends on the fugacity b. In Figure [2] 
we show for each component and each scenario the clan parameters as a function of the fugacity. 
Again, the points come from fits to experimental data, the solid lines are our extrapolations. The 
dashed grey lines show the variation of clan parameters with b at fixed k (that is, at fixed V and 
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b = fugacity = n/( n + k ) 



a = k b 
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Figure 1: Fugacity and a parameter dependence on cm. energy 
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N = avg number of clans 
= Hn(l -b) 



n c = avg num of particles per clan 
= [ (1-1/b) ln(l-b)] 




Figure 2: Clan parameters dependence on fugacity. 
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T) for the following values of k: 1 (lowest curve), 3, 7, 30 (highest curve); being n c independent 
of k, only one dashed line is visible in the corresponding graphs. 

For the soft and scenario- 1 semi-hard components it is shown that N is a very slow growing 
function of the fugacity of the system throughout the ISR region (b < 0.7), and then starts to 
grow quickly; n c as a function of the fugacity has a similar behaviour from « 1.5 to ~ 3. 

The decrease of the average number of clans in scenarios 2 and 3 leads again to the conclu- 
sion that this behaviour is closer to that of a quantum system than to a classical one, favouring 
as it does the production of larger clans and therefore of regions of higher particle density. 

For a discussion of other parameters, like the void scaling function, see again Ref. [j^]. 

4 Conclusions 

By defining a new thermodynamical approach to multiparticle production at parton level we 
have given the physical meaning of fugacity to a parameter previously used only to describe 
deviations from Poisson behaviour in multiplicity distributions. 

On this basis, we revisited our previous extrapolations to the TeV region of inelastic hadron- 
hadron collisions and examine the different behaviours of the two classes of events ('soft' and 
'semi-hard'). 

In the first class, i.e., soft events, the ancestors of the clans are produced earlier, at higher 
virtuality and when the temperature is higher. Ancestors in these conditions generate little (clans 
are small). This results in a behaviour closer to that of a classical system (ideal gas). 

In the second class, i.e., semi-hard events, the ancestors are produced later in the cascade, 
at lower virtualities and when the temperature is lower. Ancestors in these conditions are more 
prolific (clans become larger). This results in a behaviour closer to that of a quantum system 
(stimulated emission); high density regions exist. 

Although we used explicitly in the illustration the NB(Pascal)MD, our result is extensible in 
principle to any infinitely divisible distribution which also belongs to the class of power series 
distributions. 

The results discussed in this paper bring in the spotlight the concept of clans, which up 
to now was only applied in a statistical framework. At this point, it becomes important to 
investigate other physical properties of clans, in order to answer questions like the following 
ones: can clans be considered observable objects? if so, what are their quantum numbers? do 
they start to interact among themselves in the TeV region? how will this possibility modify the 
ideal gas equation of state? 

Work in this direction has already begun[Q| by studying clan masses, with preliminary indi- 
cations that the answer to the first question (observable clans) is positive. This can be extremely 
relevant for the new heavy ion machines where the standard examination of events with tens of 
thousands of particles may be very problematic. 
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